Pressure Corrections to the Equation of State in the Nuclear Mean Field. 
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We show the connection between stiffness of equation of state in a Relativistic Mean Field (RMF) 
of Nuclear Matter (NM) and the existence of a strong violation of longitudinal Momentum Sum Rule 
(MSR) in RMF for a finite pressure. The increasing pressure between nucleons starts to increase 
the ratio of nucleon Fermi to average single particle energy and according to the Hugenholtz-van 
Hove theorem valid for NM the MSR is broken. In order to satisfy that MSR we propose changes 
which modify the nucleon Parton Distribution Function (PDF) above a saturation density. The 
course of Equation o State in our modified RMF model is very close to semi-empirical estimation 
and to results obtained from extensive DBHF calculations with a Bonn A potential. Other features 
of the model includes a good values saturation properties including spin-orbit term. The proper 
stiffness of EoS recently discussed in an application to compact and neutron stars is important when 
studying star properties (mass-radius constraint), especially the mass of "PSR J16142230" the most 
massive known neutron star, which rules out many soft equations of state including exotic matter. 
An admixture of additional hyperons are discussed in our approach. 

PACS numbers: 24.85.+p 



I. INTRODUCTION 

Experimentally, in Deeply Inelastic electron Scattering 
(DIS) on nuclear targets, fotons with large negative mo- 
mentum square — q 2 = Q 2 > lGeV 2 and large energy 
transfer v, probe interacting hadrons - a kind of moving 
sub-targets. Start with the picture of a nucleus with mass 
Ma- Bjorken scaling allows to describe nuclear dynam- 
ics by the Parton Distribution Function (PDF) F^(xa) 
which depends on the Lorentz invariant Bjorken vari- 
able x A = Q 2 /(2M A v)^. Generally the SF depends 
also on the resolution Q 2 which is particularly important 
for x a < 0.01 1 A where nuclear shadowing takes place. 
Shadowing should be included in any treatment of the 
EMC effect. However shadowing is described0-01 as a 
multi-scattering process with diffraction between differ- 
ent nucleons. If the Momentum Sum Rule (MSR) has 
to be analyzed, the simple convolution of nucleon PDF 
with nuclear distribution preserves the Longitudinal Mo- 
mentum (LM) of this parton system. For xa > 0.1/ A we 
know0,|j-0] that nuclear shadowing is unimportant. 

In the light cone formulation^, Q, xa corresponds to 
the nuclear fraction of quark LM k + = k° + k 3 and 
is equal (in the nuclear rest frame) to the ratio xa = 
k + /Ma- But the composite nucleus is made of hadrons 
which are distributed with longitudinal momenta p^, 
where h = N, n, ... stands for nucleons, virtual pions, ... . 
In the convolution modelJfllH a fraction of parton LM xa 
in the nucleus is given as the product xa = Xh*Uh of frac- 
tions: parton LM in hadrons Xh = Q 2 / {2Mh,v) = k + /p^ 
and LM of hadrons in the nucleus y% = pX/Ma- The 
nuclear dynamics of given hadrons in the nucleus is de- 
scribed by the distribution function fh(y = Vh) and PDF 
FJ/(x = Xh) describes its parton structure. Let us re- 



member that there are two different scales of interaction: 
long range nuclear scale which forms hadron distribution 
functions in nuclear matter and a much shorter parton 
scale which is responsible for their PDF's. In the con- 
volution model restricted to nucleons and pions (lightest 
virtual mesons) the nuclear PDF F 2 is described by the 
formula: 

F?{xa) = JydyJdx5(xA-xy)(f N (y)F 2 B (x) + r(y)FZ(x)) 
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d 4 p A(p°+p 3 ) 



(27T) 



^i)Tr [ 7 +5(p,P)](l) 



where KT and 



F 2 are the parton distributions of the 



virtual pion and bound nucleon, P is the total four mo- 
mentum of the nucleus. The function S(p, P) is pro- 
portional to the connected part of the nuclear expecta- 
tion value of the nucleon Green's function Q and the 
trace is over the Dirac and isospin indices. Both the 
quark and nucleon distributions in the the formula ([1]) 
are manifestly covariant and can be expressed by Green's 
functions^. The single nucleon Green's function S(p, P) 
in nuclear medium is given for example in 0, H| : 



S( P ,P) = -i( 1 -(p-U v ) + M* N )x 



(2) 
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E*(p) 



M 



N 



6(p°-E*(p)-U°) 



(3) 
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The values of vector XJy = gv{Vo,0) and scalar (M^ — 
M/v) potentials are discussed in three specific mean-field 
models [71 y, [l0(. In the nucleus rest frame V~ = V + = 
V° for the expectation values of vector meson fields. The 
connected part (second term) of Q inserted into Eq. (TTJ) 
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for /iv(y) gives, after taking the trace and using the delta 
function to integrate over p°, the result which can be 
simplified in the RMF by using the Hugenholtz-van Hove 
theorem to the form [ill: 



4 f S N (p)d 3 p p 3u/ +/ . 

—3— ( 1 +-E^) s (y-p + /za) 
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Here the nucleon spectral function was taken in the im- 
pulse approximation: Sn = n(p)8(p° — (E (p) + E/y))> 
Ea = Ma/ A. Ep is the nucleon Fermi energy and y takes 
the values given by the inequality (Ep ~pp)/ea < y < 
(Ep +pf)/ea- The Bjorken scaling in xa corresponds 
in light cone dynamics to the scaling with respect to the 
product xy. This scaling assume that due to relativis- 
tic contraction we can in first approximation neglect the 
qq and AW interaction in such a deep inelastic parton 
structure. 

The MSR for the nucleonic part is sensitive to the 
Fermi energy as can be seen from the integral: 



dy yfN(y) 



El 

£A 



(5) 



Thus the nucleonic part of MSR gives a factor Ep/sA 
which is equal to 1 at the saturation point. 



F A (x)dx = / dyyf N (y) F 2 B ' (x)dx= F 2 B (x)dx (6) 



To describe the "EMC effect" at the saturation where 
— = 1 the nucleon PDF should be gradually smaller [1, 
[3, [l3| in medium for Bj0"rken x greater than the single 
nucleon resolution xl ~ 1/Af/vP 3 ~ 0.4. Consequently 
with eqQ] the average pion excess < p-n + > is given by 
the difference: 



II. NON-EQUILIBRIUM CORRECTION TO 
NUCLEAR DISTRIBUTION. 



For finite pressure very important is well known 
Hugcnholtz van Hove relation between Ep , ea and pres- 
sure p (see for example [lT}) which was proven in self- 
consistent mean field approach [l8| to NM. The Fermi 
energy is defined as density derivative of the total nu- 
clear energy E A = Aea- 



jpA _ d (jpA\ d 



£a _|_ q ^ £a (9) 
dg V Q J dg 



where A/ g = Q gives the volume. The external pressure 
Pe is given by the thermodynamical relation: 

fdE A \ yd,, 

PE ^ -(«rJ A = ^^ (10) 



E F 



£A + Et 



(11) 



Where the pressure dependent part E press = p/g. 

The condition ([5]) which measured the fraction of the 
nuclear LM carried by partons in nucleons gives for neg- 
ative pressure pe- 

f dy yf N (y) = ^ = 1 + p E /g < 1. (12) 
J ?A 

The possible interpretation: the missing part of LM 
in (|12|) is carried by quarks localized inside attractive 
mesons(pions) responsible for this negative pressure. 

However we are interested in high density region with 
positive pressure due to a nucleon repusion. Above the 
saturation point the increasing pressure between nucleons 
starts to increase the nucleon Fermi energy Ep so MSR is 
broken by the positive factor Ep/sA- It spoils the mean 
field approximation and indicates that increasing pres- 
sure distorts the parton picture where the constituents 
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F A (xA)dxA - 
f (F 2 B (x)-F 2 N (x))dx 

J XT, 



Fy (x) dx 



(7) 

(8) 



and is sufficiently small d, [ID, [3 (~ 1%) to describe 
also the nuclear Drell-Yan reactions. 

However in our interacting system part of nucleons oc- 
cupy states above the Fermi level. Therefore our formula 
(0| and MSR should be treated as the first approxima- 
tion. Alternatively the mean field scenario should be 
supplemented by neutron-proton Short Range Correla- 
tion (SRC) which has the remarkably similar A depen- 
dence as the EMC effect fl(| . On the other hand the 
good description of the nuclear Drell-Yan reaction with 
the sufficiently small [1, [2 [HI (~ 1%) pion admixture 
provides that our convolution approximation is up to few 
% accurate. 
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FIG. 1: R(x) - The density dependent energy E pre33 inside 
NM for stiff (Walecka) and soft (ZM) EoS. 



3 



s are weekly bound in the large momentum transfer can 
be treated as free. Treating nucleons as bag s the finite 
pressure will influence their surfaces [19l - l2]| . To com- 
pensate this the nucleon PDF in the nuclear medium has 
to be adjusted. In particular it will be shown how the 
nucleon mass in medium should start to decrease in re- 
sponse to positive pressure. In the paper we show how 
such a relation change the EoS. Consider the nuclear 
pion contributions above the saturation point. The av- 
erage distances between nucleons are smaller, therefore 
parameter Xi (|5|) increase with density. It means that 
the room for nuclear pions given by (J5J will be reduced 
for higher densities. Also the pion effective cross section 
is strongly reduced at high nuclear densities above the 
threshold in A + A = A + A + 7r reaction calculated 
in Dirac-Brueckner approach I22f ( also with RPA inser- 
tions to self energy of N and A [23[ included) . Therefore 
for positive pressure the nuclear pions carry substantially 
less then 1% of the nuclear LM and dealing with a non- 
equilibrium correction to the nuclear distribution (JXJ) we 
will restrict considerations to the nucleon part without 
additional virtual pions between them 

jF 2 A (x)dx = jj dyf N (y) J F 2 N (x/y) . (13) 

The eventual admixture of additional pions makes the 
violation of longitudinal momentum even stronger. The 
Equation of State (EOS) for NM has to match the sat- 
uration point with compressibility A' -1 = 9q 2 j^2^x ^ u t 
then the behavior for higher densities is different for dif- 
ferent RMF models. Generally, the choice of initial La- 
grangian in the nuclear RMF models and the depen- 
dence of nucleon masses from density is not unique. 
We compare here two extreme examples: stiff model 
of Walecka Ml a nd soft non-linear Zimanyi-Moszkowski 
(ZM) model 111 H3] (sea also [H[). For linear coupling 
in the standard Walecka model at saturation density of 
NM compressibility is too large (if -1 ~ 560 McV). The 
Moszkowski-Zimany [Toj model produces the very soft 
EoS (see marked dotted line in fig|2]) with a good value 
of A'— I = 225MeV. The energy E press = pe/q shown in 
FigJT] calculated for those two models reflects courses of 
the nuclear EOS. We see from that the E press is relatively 
small in the ZM model (for g = 0.3fm^ 3 only 25 MeV 
in comparison to 150 MeV in Walecka model). In both 
models, two coupling constants of the theory are fixed 
by the empirical saturation density. The additional (to 
the Walecka model) derivative coupling of the nucleon to 
the scalar field enables to re-scale the fermion wave func- 
tion, and interprets the new, density dependent nucleon 
"Dirac" mass M z m 0, El HI . It starts to decrease with 
density, and at the saturation point reaches 85% of free 
nucleon mass. But from the point of view of our descrip- 
tion of the EMC effect [J El El, this means that 15% 
of nucleon longitudinal momentum is missing, and there- 
fore has to be compensated by the enhanced pion cloud. 
Such a big pion cloud, however, spoils the description of 
nuclear Drell-Yan experiments (a. ff2l. [l5| which measure 



the sea quark enhancement. There is no evidence for a 
such huge enhancement in the EMC effect for x = 0.2. 
At the density p = 3p = 0.48/to" 3 M Z m ^ 0.73Mjy. 

Our approach has different attitude, the nucleon mass 
begin to decrease above the saturation point. The grow- 
ing pressure between nucleons starts to increase the Ep 
([5"j) and consequently the sum rule is broken by the 
factor Ef/sa > 1- To compensate this factor which in- 
creases the longitudinal momentum (|6]) of nuclear par- 
tons, the the nuclear distribution and the nucleon PDF 
in the nuclear medium has to be changed. For good es- 
timate, in order to proceed without new parameters, as- 
sume that the changes of PDF will be included through 
the changes of nucleon mass Mjy in NM like in a QCD 
motivated bag model(20|. Multiplying the argument of 
the PDF by a factor Ef/sa the PDF will be squeezed 
towards smaller x and the total fraction of longitudinal 
momentum will be smaller by a factor e^/A^: 

[ F 2 N f^a*) dx N = £A[ Tf F 2 N (x)dx^^ [ F 2 N (x)dx 

JO \ £ A J &FJo E f J 

(14) 

Here in the integral we neglect the small contributions 
from x > 1 region originated from NN correlations. 

Now, with the help of Eq.© and Eq . (fl"4"|) . the nuclear 
MSR is satisfied: 

-j jF 2 A (x A )dx A = j dyyf N {y) J F 2 \ J^- d XN ^j dxN 

This means that quarks in the nucleus carry the same 
fraction of longitudinal momentum as in bare nucleons. 

E A 

Such a scaling of Bjorken ratio a; at by the factor in 
the integral (fT4|) changes the total sum of the quark lon- 
gitudinal momenta k~^ i = k a Ni + kf !i proportional to the 
quark energy in the nucleon ( nucleon mass Mn). Con- 
sequently the nucleon mass in medium M me d for g > go 
will decrease with the density by the gradually decreasing 
factor SA/Ep-. 

M med =^M N =M N / ( l + 9 AM) ^M N (l-PE.) 
Ep \ e A J \ ge A ) 

(15) 

The generalized Hugenholz van Hove theorem can be 
obtain [27J for different kinds of baryons in matter, for 
example S strange barions (E, A) in a nuclear matter (28j 
with total barion number - (A+S). Analogously to Eg ITU 
a sum of all Fermi energies, including Fermi energies of 
strange barions Ep, is equal to the total energy 

E A+S = 

Ma+s with analogous pressure corrections: 

AE^ + SE^ = Ma+s + {A+S) PE /g (16) 

fdE A+s \ 
with p E = - — — — 

V dn J A+S 



4 



The corrections (|15l) depend only from pressure and total 
energy of the system therefore the similar medium cor- 
rection based on Eq[l6]will apply to masses of different 
barions, like strange A and S, in NM. 

Let us discuss these results in the simple bag model [29j 
where the nucleon in the lowest state of n = 3 quarks is 
a sphere and its energy E Ba g is given as a function of 
radius R with phenomcnological constants - w , Z and 
B: 



EBaq 



nuj ~Z t 4tt 3 



R 



(17) 



This formula was obtained by studying the free quark 
flow described by the stress tensor Tp" is given on bag 
surface along the vertical to this surface by: 



n - 



l_d_ 

2dx h 



1a( x )Qa(x) 



(18) 



where q a (x) denotes the quark field with given flavor and 
color. Because the q a (x)q a (x) = on the surface its 
derivative must lie along the normal and define the in- 
trinsic pressure po- 



n„PD 



(19) 



The pressure condition for can be obtained from a 
usual critcrium (|11[) : 



{dE Bag /dn Bag ) n = o 



(20) 



However in a compressed medium the pressure pc gener- 
ated by free quarks inside the bag is balanced at the bag 
surface [29] not only by a intrinsic confining force repre- 
sented by the bag constant B but additionaly by exter- 
nal pressure pe generated by elastic collisions with other 
hadrons (l9l [20L |30| . In this way the residual strong inter- 
action with all other quarks localized outside the given 
bag is taken into account : 



p D = B+pe 



(21) 



Finally using the equation Eq. (|20p and resulting spher- 
ical solution for the bag radius we obtain: 



3lUq — Zq 



R 2 



R = 



A-n(B+p E )R 2 

3u>q — Zq 



"I 1/4 



4?r(B +pe) 



(22) 



At the saturation p B = and the bag constant B 
can be determined by the value of the nucleon ra- 
dius R ~ 0.6. fro. Above the saturation point the R 
increases (20l [2lj] . The main reason: the external pres- 
sure was not taken into account and the value of B 
wrongly determines the increasing bag radius R (nucleon 
swelling). Here we include explicite the external pres- 
sure which increases when the bag constant decreases. 
When the sum (B + pe) weakly depends from density, 



then the bag radius will remain, according to (|22|) ap- 
proximately constant with the increasing nuclear density 
q. In fact such a assumption is confirm in the calcula- 
tion of [20| where an decreasing of the B constant from a 
saturation density p up to 3p is accompanied be similar 
increase of the external pressure pe- The QMC model 
in medium [U [3l[ allows to find how the nucleon radius 
is changing with density and for the derivative coupling 
model with ZM[l0| Lagrangian the nucleon radius re- 
mains almost constant [3 1| up to density g = 10po- 

If for g > go B(g) + pe(q)) — B(g = go) that the 
nucleon energy in the rest frame can be written as: 



E 



Bag 



AttR 3 



E 



PE 



Bag 



(23) 



with g ~ 1/Q.Bag > go where the volume of a space 
between nucleons above the saturation density go was 
neglected. E Bag differs from the the nucleon mass by 
the cm. correction [2l| to the nucleon parton model. 
Now we can compare the expressions (fT5)) with (f2Uf for 
the nucleon mass in medium. The pressure corrections 
essentially given by E press = (p/g) are similar, espe- 
cially when the nucleon radius and (B + p B ) remains 
constant 3l| in the bag (|2"3"|) model. Such a solution of 
R(g) which is constant with g testifies the properly soft 
EoS; i.e. ZM model with K' 1 ~ 250 McV showed 
in Figs. 1,2. The pressure corrections to the masses of 
strange baryons will be different because the important 
ingredient the mass of the strange quark( 300MeV) re- 
mains constant and the reduction ([23]) concerns the re- 
maining masslcss part. In effect the mass difference be- 
tween strange baryons and nucleons will increase slightly. 
Also the best candidate ]^ . |32| to strange admixture 
the E~ has the rcpulsive[33] S — NM interaction which 
also shifts the possible strangeness admixture to higher 
density [32j]. Therefore we present now the calculation re- 
stricted to nucleons as a good estimate for densities from 
go to 3 go. The full calculation which will include admix- 
ture of strangeness will be published elsewhere. 

Now our explicit mass dependence from density, en- 
ergy ea and energy derivative (fT5|) is plugged into the 
following standard Walecka RMF equations [7| for energy 
per nucleon ea and effective mass M*\ 



£A =Clp+^(M med -M* med ) 2 +—l— [d 3 pJ(p 2 +M* 2 ) 



Ml 



2C 2 2 (27r)3 



d 3 p 



M* 



\I{P 2 + M* 2 ) 



(24) 



where 7 denotes level degeneracy (7 = 2 for neutron 
matter) and two (coupling) constants: vector C 2 and 
scalar C 2 , were fitted p], |9( at the saturation point of 
nuclear matter (in the formula 2C 2 — C 2 /M j 



2C 2 = M 2 N /C 2 ). 



and 
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FIG. 2: The nucleon energy ea — Mn as a function of NM 
density for two RMF models; scalar-vector Walecka (dot lines) 
and our Modified Mass approach (solid). Both RMF models 
are calculated for two parameterizations: Si version[3| (po = 
.19/m-f) and version[<| S 2 {po = .16/m -3 ). Results for full 
DBHF[36(] (dotted marked line) calculation using Bonn A NN 
interaction are displayed for comparison, also nucleon energy 
m ZM modelQI} is in the plot (dotted marked line). 
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Note that equation (|24[) is obtained from the energy- 
momentum tensor for the model Hamiltonian with con- 
stant nucleon mass @ ■ Here we assume that the same 
equation with medium mass M mec i will be satisfied. It 
should be a good approximation, at least not very far 
from the saturation density. The pressure p is obtained 
from the thermodynamic relation 

The final results were obtained by solving 1 numeri- 
cally differential recursive equations (|25p. starting from 
standard solutions of Eq.(|24p at the saturation density 
for two versions of the Walecka model: first version 
Sipf have a minimum at po = .19/m~ 3 (parameters 



0.2 0.3 0.4 0.5 0.6 



C v z = 195.9, C/ = 267.1) and second version S 2 @ 
(po = .16/m^ 3 ,C„ 2 = 273.8, C s 2 = 357.4). They arc 
displayed in Figs.([5]|3|. In FigJ2]our values of energy per 
nucleon ea calculated for two versions S\ and S2 are de- 
noted by solid lines (Mass Mod.) and solutions of the 
ordinary Walecka model with constant mass, denoted by 
dashed lines, are presented for comparison. Our EoS's 
are generally much softer - from the unrealistic value of 
A' -1 = 560MeV for the Walecka model (Si) to the rea- 
sonable K— 1 = 290MeV obtained in our model. Below 
saturation density the solutions are of course identical 
(solid lines). The nucleon mass M me d and the nucleon 
effective mass M* (both in Walecka and in our (Mass 
Mod.) model) are compare in FigJ3] The decrease of 
M me d is directly connected with the change of nucleon 
SF for positive pressure in NM. The decrease of M me d 
above saturation is directly connected with the change 
of nucleon PDF for positive pressure in NM. In FigH] 
we plot the pressure against NM density. Our results 
are compared here with a semi-experimental estimate [35[ 
from heavy ion collisions and indeed correct (solid lines) 
Walecka results (dashed) quite well making the EoS sig- 
nificantly softer. For SI parametrization our model in 
FigHJfits well the allowed region in neutron matter. For 
nuclear matter we have good course of EOS below den- 
sity p = 5/m~ 3 . In fact, for this density, the (partial) 
de-confincmcnt is expect which will change EOS above 
the phase transition [33]. Our results are close (slightly 
below) DBHF results (dotted line). Therefore it is in- 
teresting how the realistic NN calculations with off-shell 
effects correspond with the modification of the nucleon 
PDF in medium induced by MSR on the parton level. 



FIG. 3: Constant nucleon mass Mn (Walecka) and density 
dependent mass M me£ j from our "Mass Mod." model. Also 
effective mass M* in these approaches. Both models are cal- 
culated for Si and S2 parametrization. 

In the Walecka model M me d = Mn- In our model the 
finite pressure corrections to Mmed (fT5j) convert the re- 
cursive equation (|2"4")l to a differential equation above the 
saturation density po in general form: 

f(e A , (e A )) = Oforp> Po (25) 



It is important to mention that in solutions of Eas J24l25ll the 
Fermi energy from definition Eq. J9} has a different value from 

the usual form E F = ^ M* 2 + p% + U v used in Eq.Q. The 
discrepancy vanish near the saturation density, increases with 
the density and reach the 15% of the total vector repulsion 
in EqJ24t . Similar problems are well known in the saturated 
NM [3411 and are connected with the proper choice of single parti- 
cle potential which in our case should be adjusted to the changes 
of nucleon mass. This discrepancy can be removed, here e.q. by 
the less repulsive momentum dependent vector potential for nu- 
cleons however such a correction has no influence on presented 
mass changes of the nucleon and our final results. 
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For example, in DBHF method, there are additional 
corrections [36j from self energy which diminish the nu- 
clcon mass with density. 

III. CONCLUSION 

We have presented a simple, parameter free, hybrid 
model of nuclear matter, which examines the influence of 
the parton structure on nuclear EoS. It was shown that 
there is connection between stiff EOS and strong viola- 




p (fm-3) 




FIG. 4: Pressure for neutron matter (upper plot) and for nu- 
clear matter (bottom) as a function of density for two most 
frequent pure scalar-vector parameterizations of Walecka 
RMF model: version Si {po = .16/m" 3 ) and version^ 
S2 (po = .16/m -3 ) are denoted by dashed lines. Our (Mass 
Mod.) modification are denoted by solid lines. The area 
denoted by "flow constrain" taken from [35| determined the 
allowed course of EoS, using analysis which extracted from 
matter flow in heavy ion collisions the high pressure obtained 
there. The DBHF fref. [36(|) calculation with Bonn A interac- 
tion are shown as dotted line for comparison. 



tion of longitudinal momentum conservation carried by 
quarks. In our approach we propose a scaling of the nu- 
cleon mass which modify the nucleon PDF in order to 
satisfy the longitudinal momentum sum rule for the par- 
tons in NM above the saturation density. Our pressure 
correction reduce the violation of MSR from 50% in the 
linear Walecka model to 10% in our model. This modify 
the EoS for nuclear and neutron matter making it softer, 
close to semi-empirical analysis (HJ and close to DBHF 
calculation with a realistic Bonn A potential [3(|. The 
similar mass corrections concern all baryons in medium 
and will be present e.q. in the strange nuclear matter. 
Other features of the Walecka model, including a good 
value of the spin-orbit force remain in our model un- 
changed. Our results suggest pressure corrections above 
the saturation density to any RMF model with constant 
nucleon mass and unmodified parton PDF (without a 
properly normalized nucleon longitudinal distribution). 
Finally we argue that the nucleon radius or the con- 
fining region should remains constant with density oth- 
erwise the pressure correction to the nucleon mass will 
be different in Bag model Eq. ([2"3")) versus parton model 
Eq. (ll5[) . The alternative: increasing confining region in 
the Bag model which corresponds eventually to the in- 
creasing number of intermediate pions(mesons) for a high 
density will give a very stiff EoS; rule out when studying 
star properties (mass-radius constraint), especially the 
mass of PSR J16142230 the most massive known neu- 
tron star |3^| recently discussed[37| in the application to 
compact and neutron stars. Partial support of the Min- 
istry of Science and Higher Education under the Research 
Project No. N N202046237 for is acknowledged. 
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